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We study the statistics of charge transport in a mesoscopic three-terminal device with one super- 
conducting terminal and two normal-metal terminals. We calculate the full distribution of transmit- 
ted charges into the two symmetrically biased normal terminals. In a wide parameter range, we find 
large positive crosscorrelations between the currents in the two normal arms. We also determine the 
third cumulant that provides additional information on the statistics not contained in the current 
noise. 



The number of charges transfered in a transport pro- 
cess fluctuates due to quantum-mechanical uncertainty 
and statistics. Therefore, the outcome of a current mea- 
surement accumulated over some time period to is in gen- 
eral described by a probability P{N), where N is the 
total number of charges transfered. P{N) is called the 
full counting statistics (FCS) of the transport process [Q . 
The first two moments of the FCS are related to the av- 
erage current and the current noise and are accessible 
to present experimental techniques. Higher-order corre- 
lations are likely to be measured in the future. Several 
schemes to measure either higher correlators or the full 
distribution have been proposed recently ||, |[ 0, ^, || . 

The current noise, i. e., the second moment of the FCS, 
is of particular interest. It can be used as a diagnostic 
tool to probe the nature and the quantum statistics of 
the charge carriers and the existence of entanglement 
1^. For superconductor(S)-normal metal(N) heterostruc- 
tures, a doubling of the shot noise in comparison to the 
normal case was predicted and measured in diffusive 
heterostructures . Recent calculations taking into ac- 
count the proximity effect in such structure s jllj are in 
good agreement with experimental results |12|]. Multi- 
terminal SN structures have been suggested to produce 
entangled electron pairs Q . 

So far, crosscorrelations, i. e., current correlations in- 
volving different terminals, were measured only in nor- 
mal single-channel heterostructures . These have con- 
firmed the prediction that current crosscorrelations 
in a fermionic system are always negative. To our knowl- 
edge, there is no measurement of crosscorrelations in a 
system with superconducting contacts up to now. The- 
oretically, positive crosscorrelation with a single-channel 
beam splitter for Andreev pairs injected from a super- 
conductor have been predicted ||l7|| . In a setup in which 
crosscorrelations between a normal lead and a tunnel- 
ing probe are considered, the sign of the correlations was 
found to depend crucially on the sample geometry [ p^ . 
A numerical study found positive crosscorrelations in a 
three-terminal device with a few channels with ferromag- 
netic contacts 

In this Letter we find the full counting statistics of a 
many-channel beam splitter that divides a supercurrent 
in two normal quasiparticle currents. We calculate the 



distribution of the transmitted charges taking the prox- 
imity effect into account. For comparison we also calcu- 
late the FCS for the case in which the superconducting 
terminal is replaced by a normal one. 

The setup of our three-terminal device with one su- 
perconducting and two normal-metal terminals is shown 
in Fig. |l|. All three terminals are connected by tunnel 
junctions to a small normal-metal island. We assume the 
island to be large enough, that we can neglect charging 
effects, and small enough that we can treat the Green's 
functions of the island as constant. Thus, we are re- 
stricted to energies below the Thouless energy of the is- 
land. The system is then appropriately described by the 
circuit theory of mesoscopic transport |20| , within which 
the counting statistics is easily accessible [|[ ^ ^ . 

The circuit elements that are important for our case 
are normal, resp. superconducting terminals and tunnel 
connectors. The terminals are described by quasiclassi- 
cal Green's functions, which depend on the type of the 
terminal (N or S), chemical potential, temperature T and 
a counting field %. We assume zero temperature and a 
symmetric bias at the two normal terminals. Since we 
are interested in correlations between currents in differ- 
ent terminals, we introduce different counting fields. The 
voltage is chosen such that eV <C A. Charge transport 
(at r = 0) occurs then only in the interval \E\ < eV and 
we need to consider only this energy interval below. 

The Green's functions of the two normal terminals are 




FIG. 1: Three-terminal beam-splitter. One superconducting 
or normal terminal (S or N) and two normal terminals (A'^i 
and N2) are tunnel-coupled by conductances g, gi and (?2 to 
a common central node . A current is passed from S/N into 
the two normal leads that are kept at the same voltage. Ideal 
passive charge counters are indicated by the counting fields 
X, Xi and X2. 
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then given by 

Gl,2 = e''^^arK|1Q^^-^xu^ft,|2 



(1) 



where Gat is the same for both normal terminals. At 
zero temperatm^e Gjv = ctsts + (fi + if2)l for \E\ < 
eV and Gat = 0-37=3 + sgn(i?)(T3(fi + 1x2) for \E\ > eV. 
Here (Ti(fi) denote Pauli matrices in Nambu(Keldysh)- 
space. The comiting rotation matrix is tk = ^3Ti. The 
superconducting terminal in equilibrium is characterized 
by Gg — &1I and a counting field Xi that enters as in (|^). 

If one node is connected to M terminals by means of 
tunnel connectors, one can find a general form of the 
FCS, i.e. the probability P{Ni, ... , Nm) that iVi(2,... ,m) 
charges are counted in terminal 1(2,... ,M). The un- 
known Green's function of the central node is denoted by 
Gc. The matrix currents into the central node are given 
by /fc = ^ [Gc, Gk] , where the index k = 1, . . . ,M labels 
the terminals and gk is the conductance of the respective 
junction. The Green's function of the central node is de- 
termined by matrix current conservation on the central 



node, reading YlkLi h 



YlkLl 9kGk, Gc 



0. Em- 



ploying the normalization condition G^ ~ 1, the solution 



Gc 



Efc^m=l 9k9m {Gk, G,„} /2 



(2) 



To find the cumulant-generating function (CGF) 
S of P{Ni,...,Nm) we integrate the equations 
{~ito/e)dS{xi, XM)/dxk - / dETiTKlk/Se [||. 
We obtain 



S{xi 



to fdE^^ 
e J 2 
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{ Gk , G m } . 

(3) 



This is the general result for an M-terminal geometry 
in which all terminals are tunnel- coupled to a common 
node. 

We now evaluate Eq. for our three terminal setup. 
Introducing pi — 2ggi/{g'^ + (ffi + 52)^) we find 



r 



(4) 



This result for the cumulant-generating function incor- 
porates all statistical transport properties for our present 
setup. The inner argument contains counting factors for 
the different possible processes. A term exp(i(xfe + Xi ^ 
2x) — 1) corresponds to an event in which two charges 
leave the superconducting terminal and one charge is 
counted in terminal k and one charge in terminal The 
prefactors are related to the corresponding probabilities. 
For instance, pi is proportional to the probability of a 
coherent tunneling event of an electron from the super- 
conductor into terminal 1. A coherent pair-tunneling 
process is therefore weighted with p\. This is accompa- 
nied by counting factors which describe either the tun- 
neling of two electrons into terminal 1(2) [counting factor 
exp(«2(x]^(2) ~ x)) ~ 1] or tunneling into different termi- 
nals [counting factor exp(i(xi+X2 — 2x)) — l]. The double 
square-root function shows that these different processes 
are non-separable. 

It is interesting to compare Eq. (^ with the case in 
which the superconductor is replaced by a normal metal. 
The resulting CGF is 



^'^(xi,X2,x) = -^(.g + .gi + .g2)x 

^1 + (e«(xi-x) - 1) (e»(x2-x) - 1) , 



(5) 



where ^^[2) = ■^99i{2}/{9 + 9i + 92f- Thus, one of the 
square roots in Eq. (|^) can be attributed to the multi- 
ple tunnel-junction geometry, which is already present in 
the normal configuration. The second square root in the 
CGF for the superconducting case must then be due to 
the proximity effect. 

We now evaluate some average transport proper- 
ties of the S|NN-system and compare them to the 
N|NN-case. The currents into the different termi- 
nals are obtained from derivatives of the CGF: Ik ~ 

{—ie/tQ)dS/dxk\xi=x2=x=o- ^^'^ trans-conductances 
Gk = Ik/V into terminal fc 1, 2) are then given by 



Gk — 



^5fc(5i +52) 



(.9' + (51 



92) 



3/2 
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99k 



9 + 91+92 



■•(6) 



The superscript S{N) denotes the S|NN(N|NN)- 
case. Noise and crosscorrelations are obtained 
from second derivatives of the CGF, i. e., P^^ — 
i2e^/to)d^S{xi,X2,x)/dXkdxi\xi=X2=x=a- We define 
Fano factors Fki = P^;/2e/, and we denote the Fano 
factor of the total current with F — Fu + F22 + 
2Fi2. We also calculate the third cumulant of the to- 
tal charge transfer (normalized to the Poisson value) 
G3 = (ie//io)93S'(0,0,x)/ax^lx=o- The results in the 
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FIG. 2: Conductance, Fano factors, crosscorrelations and 
third cumulant of the beam spUtter. The thick hnes corre- 
spond to the S|NN-case and the thin hnes to the N|NN-case. 
The conductance (upper-left panel) in the superconducting 
case shows a maximum around g = g\ + 92- In the normal 
state, the conductance varies between g and (?i -f (?2- In both 
cases, the current noise (upper-right panel) shows a suppres- 
sion around g = gi + g2 as compared to the limiting values 
of 2 in the S|NN-case and 1 in the N|NN-case. Large positive 
crosscorrelations occur in the superconducting case (lower- 
left panel), whereas they are always negative in the normal 
case. Around g — gi + g2, the superconducting crosscorre- 
lations become negative. Note, that what is plotted here is 
fi2'^ ~ P'i2^ {g^ + 52)^/<?ig2. The third cumulants (lower- 
right panel) are always positive. Around g ~ gi + g2 they 
are strongly suppressed. In the S|NN-case, Cf has a double- 
minimum here, as shown in the blow-up. 



superconducting case are 



pS 9192 
^12 - , „ ^2' 



(ffi + 92^ 



(l-5a;^) , F'' =2-5x^ 



(7) 



_ A Qn^2 , «Q^4 „ _ 3(51 + 92) 



9^ + {91 + 92Y 
In the N|NN case, on the other hand, we find 



N 



9192 



(51 + 92y 



XN , = 1 - 2xN , (8) 



= 1 — 6a; jv -f- 3x% , x n 



9(91+92) 
(.9 + 91 + .92)" 



All other Fano factors can be deduced from F12 and 
F using the relations ^^.Fki = and Fki = Fik- The 
transport properties are summarized in Fig. 0. In the 



figure the crosscorrelations are plotted as /12 = f 12 (ffi + 
92)^/9192- Most remarkably, the crosscorrelations Ff2 
are positive if x is small, whereas Fl^ is always negative in 
the normal state. Here, the Fano factor F^{F^) is close 
to 2(1). Going to the regime g (91+92) suppresses 
the Fano factor f below 2(1) and leads to negative 
crosscorrelations Fj^. In the limiting case 9 — (gi + 92) 
the Fano factors are F^ = 3/4 and F^ = 1/2, and the 
crosscorrelations are — F^ = — 9192/^(91 + 92)^- 
The third cumulant is always positive, but shows a 
strong suppression around the resonant conductance ra- 
tio 9 = (gi + 92)- In the limit of small x (xn) the third 
cumulant is 4 (1), corresponding to the effective charge 
squared transfered in a tunneling process However, 
the variation with 9/(91 +92) in the S|NN-case is more 
pronounced than in the N|NN-case. 

As an interesting side remark we point out that = 
and F^ = 1 for a;^ ~ 1/5. This looks like a signature 
of uncorrelated charge transfer in units of e. However, 
the third cumulant C§ — 13/25 differs from the corre- 
sponding value for uncorrelated le-charge transfer, viz., 
C3 = 1. Thus, higher correlations show that the charge 
transfer is still correlated. 

We briefly discuss the influence of an asymmetry gi ^ 
92 of the beam splitter. The crosscorrelations are re- 
duced, both in the S|NN and in the N|NN case. How- 
ever, the positive crosscorrelations in the superconduct- 
ing state persist for all values of the asymmetry. Cumu- 
lants of the total charge transfer like the conductance, 
F^'^ and C^'^ are independent of this asymmetry. 

Using the CGF from Eq. (^), we can identify the phys- 
ical processes leading to our previous results. We have 
seen from (^ that positive crosscorrelations are found if 
9/(91 + .92) is not close to 1. Then, pi,2 *C 1 and we can 
expand Eq. in pi.2. Dropping the trivial dependence 
on X, the CGF can be written as 



'S'(xi,X2) 



toV 



e (ff2 + (.gi+.g2)2)-V2 



(9) 



The CGF is composed of three different terms, corre- 
sponding to a charge transfer of 2e either into terminal 1 
or terminal 2 (the flrst two terms in the bracket) or sepa- 
rate charge transfer into terminals 1 and 2. According to 
the general principles of statistics, sums of CGFs of inde- 
pendent statistical processes are additive. Therefore, the 
CGF (^ is a sum of CGFs of independent Poissonian 
processes. Crosscorrelations are obtained from deriva- 
tives with respect to xi and X2- Thus, the first two 
terms in (|^) corresponding to two-particle tunneling ei- 
ther into terminal 1 or 2 do not contribute. It is only the 
last term which yields crosscorrelations, and those are 
positive. Poissonian statistics are the statistics of un- 
correlated events, which in our case means all tunneling 
events are independent. Thus, a two-particle tunneling 
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event into one of the normal terminals is not correlated 
with other tunneling events and does not contribute to 
crosscorrelations, but only to the autocorrelations. The 
two-particle tunneling into different terminals, however, 
is automatically positively crosscorrelated. The crosscor- 
relations are therefore positive. 

The total probability distribution P{Ni,N2) corre- 
sponding to (|^) can be found. It vanishes for odd values 
of (A^i -|- A^2) and for even values it is 



F(iVi,iV2) 



Ni 



^2 \ rpNi rpN2 

-'-1 ^2 



(10) 



Here we have defined the average number of trans- 
fered electrons N = t^G^V/e and the probabilities 
2^1(2) = 5i(2)/(5i + 52) that one electron leaves the is- 
land into terminal 1(2). If one would not distinguish 
electrons in terminals 1 and 2, the charge counting dis- 
tribution can be obtained from (^) by setting xi ~ 
X2 — X and performing the integration. This leads 
to Pt^t(^) = exp(-7V/2)(7V/2)^/V(7V/2)!, which corre- 
sponds to a Poisson distribution of an uncorrelated trans- 
fer of electron pairs. The full distribution ^ IS given 
by Pf^^^{Ni + N2), multiplied with a partitioning factor, 
which corresponds to the number of ways how A^i -I- N2 
identical electrons can be distributed among the termi- 
nals 1 and 2, with respective probabilities Ti and T2. 
Note, that T1+T2 — 1, since the electrons have no other 
possibility to leave the island. 

In contrast to that, we obtain in the normal case for 
t^2 1 the probability distribution: 



P''{Ni,N2)^e 



-.-N.N, 



-N2 ^^2 



(11) 



iVi! iVa! ' 

Here we have abbreviated the average number transfered 
into terminal i by A^^. Thus, the distribution in the nor- 
mal case is the product of two Poisson distributions of 
charge transfers into the two terminals. In the supercon- 
ducting case such a factorisation is not possible. 

In conclusion, we have studied the full counting statis- 
tics of a three-terminal device with one superconducting 
and two normal leads. The system is biased such that a 
supercurrent is passed from the superconductor into the 
two normal leads, with no net current between the nor- 
mal leads. Thus, the device acts as a sort of beam split- 
ter. We have calculated the full distribution of transmit- 
ted charges using the extended Keldysh-Green's function 
method fully accounting for the proximity effect. Our 
main finding are large positive crosscorrelations of the 
currents in the two normal terminals in a wide parameter 
range. These should be easily accessible experimentally. 
These positive correlations originate from independent 
Poisson processes of coherent tunneling of charges into 
the different terminals. These dominate the crosscorre- 
lations, since two-particle tunneling into the same lead 



does not contribute to the crosscorrelations. We have 
also calculated the third cumulant which provides addi- 
tional information on the current statistics not contained 
in the current noise. 
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